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ABSTRACT 

In the last few decades many efforts have been made to understand the effect of spiral 
arms on the gas and stellar dynamics in the Milky Way disc. One of the fundamental 
parameters of the spiral structure is its angular velocity, or pattern speed Op, which 
determines the location of resonances in the disc and the spirals’ radial extent. The 
most direct method for estimating the pattern speed relies on backward integration 
techniques, trying to locate the stellar birthplace of open clusters. Here we propose 
a new method based on the interaction between the spiral arms and the stars in the 
disc. Using a sample of around 500 open clusters from the New Catalogue of Optically 
Visible Open Clusters and Candidates, and a sample of 500 giant stars observed by 
APOGEE, we find ftp = 23.0 ± 0.5 km s“^ kpc“^, for a local standard of rest rotation 
Vb = 220 km s“^ and solar radius Rq = 8.0 kpc. Exploring a range in Vq and Rq 
within the acceptable values, 200-240 km s“^ and 7.5-8.5 kpc, respectively, results 
only in a small change in our estimate of ftp, that is within the error. Our result is 
in close agreement with a number of studies which suggest values in the range 20-25 
km s“^ kpc“^. An advantage of our method is that we do not need knowledge of the 
stellar age, unlike in the case of the birthplace method, which allows us to use data 
from large Galactic surveys. The precision of our method will be improved once larger 
samples of disk stars with spectroscopic information will become available thanks to 
future surveys such as 4MOST. 
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The Milky Way (MW) has long been known to posses spiral 
structure, but its fundamental nature is still under debate 
today. In the classical spiral structure theory of Lin & Shu 
(1964), as well as other models which consider that the spi¬ 
ral arms are caused by the crowding of stellar orbits (e.g. 
Kalnajs 1973; Contopoulos & Grosbol 1986; Pichardo et al. 
2003; Junqueira et al. 2013), the spiral arm pattern rotates 
like a rigid body with a well defined angular velocity ftp. In 
these models Qp is usually treated as a free parameter to be 
determined by observations. However, its value has a crucial 
importance on the understanding of Galactic dynamics and 
evolution, since it determines the place of resonances in the 
disc, for a given rotation curve. In addition to this challenge, 
there are also theories of spiral arms claiming that the pat- 


1 INTRODUCTION 


tern speed and pitch angle are variable (Toomre 1981), that 
spirals are transient phenomena on a rotational time scale 
(e.g., Sellwood & Binney 2002), or that there exist at any 
time several spiral sets with distinct pattern speeds overlap¬ 
ping in radius e.g., Masset & Tagger 1997; Merrifield et al. 
2006; Quillen et al. 2011; Minchev et al. 2012, or even that 
the spiral arms are stochastic phenomena e.g., Patsis 2006. 


* E-mail: tjunqueira@aip.de 


Looking at the existing theories, we can see that it is 
not clear whether the ftp can be described by a multiple, 
transient or a constant pattern speed. However, some au¬ 
thors have shown that the corotation radius is close to the 
solar orbit (Marochnik et al. 1972; Creze & Mennessier 1973; 
Mishurov & Zenina 1999; Dias & Lepine 2005). Amores et 
al. (2009) associated a gap in the Galactic HI distribution 
close to the corotation radius, while Scarano & Lepine (2013) 
showed a break in the radial metallicity gradient close to the 
corotation radius for many external galaxies. This suggests 
that we have dominant spiral arms with a constant pattern 
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speed for, at least, a few billion years, which do not support 
models with transient spiral arms that survive for only a few 
galactic rotations at the solar radius. Thus, the approxima¬ 
tion of a single and steady pattern speed is a useful first step 
to see whether this assumption is consistent with available 
data. 

The most direct method to measure the pattern speed 
of the MW relies on the birthplaces of the observed open 
clusters. It is done by integrating backwards in time their 
orbits according to their known location, ages and circular 
velocity in the disk. Assuming that the open clusters are 
born in spiral arms, the distribution of birthplace for some 
age bins should be spiral-like (where we assume logarithmic 
spirals), and by comparing the spiral patterns obtained from 
different ages bins, the rotation rate of the spiral arms can 
be estimated. This is a valuable way to measure the pattern 
speed, but determining ages for a large open cluster sample 
is not an easy task and may introduce large errors. 

In the present work, we introduce a new method to de¬ 
termine the angular velocity Qp based on the interaction be¬ 
tween the stars and the spiral arms, that cause an exchanged 
of energy and angular momentum (we use energy and an¬ 
gular momentum per mass unit, but for simplicity we keep 
the terminology energy and angular momentum). The only 
assumption made in this method is that the initial energy 
is the energy of a circular orbit placed at the mean radius 
Vm ■ We do not give any information a priori about the spiral 
arms, which is already included in the velocities components 
of each object. Moreover, we do not need to know the ages 
or make any assumption about the spiral arms shape. In 
addition, the fact that we do not have to make use of ages 
allows us to use observational data that provide only posi¬ 
tions and kinematics information, and hence can be a pow¬ 
erful method in the era of large spectroscopic surveys such 
APOGEE-2 (as part of SDSS-IV - Sobeck et al. (2014)), and 
future very-large ones as 4MOST^ and WEAVE^. 

The organization of this paper is as follows: in Section 2 
we present the new method for measuring the spiral pattern 
speed of the Galaxy. In this same Section we provide the 
rotation curve used to compute the Galactic potential and 
how the mean radius is calculated for each object, as well 
as some tests made to validate the method. In Section 3, 
we describe the data used to compute Qp. In Section 4, we 
compare our result with the ones found in the literature and 
discuss possible sources of errors. Concluding remarks can 
be found in Section 5. 


2 METHOD 

As a star travels around the Galactic center it interacts with 
the spiral arms, which results in exchange of angular momen¬ 
tum and energy. The rate at which this happens depends on 
the relative angular velocity Q — Qp, where Q is the observed 
stellar angular velocity. Thus, in an inertial frame of refer¬ 
ence, the energy of a star varies due to perturbations caused 
by the presence of the spiral arms but, if there is only one 
pattern speed, we can find a Hamiltonian that is time in¬ 
dependent. This Hamiltonian system lies on the frame of 

^ 4-m multi-object spectrograph telescope - de Jong et al. (2012) 
^ http://www.ing.iac.es/weave/index.html 


reference of the spiral arms, which is a non-inertial frame 
of reference. This is known as the Jacobian of the system 
(Binney & Tremaine 2008), written as: 

H = E- QpJ, (1) 

where H is the energy in the non-inertial frame of reference, 
while E and J are the energy and the angular momentum 
in the inertial frame of reference, respectively. The energy 
of star E is given by: 

+ +‘3E>o(r) + <hi(r,0), (2) 

and 

J = r'^Q, (3) 

U is the radial velocity toward the Galactic center, r is the 
star’s distance from the Galactic center, Q = V/r where V is 
the rotational velocity, 4>o(r) is the axisymmetric potential 
and 'I>i(r, 0) is the perturbative potential due to spiral arms. 
Once the total energy H is conserved {AH — 0) we can 
derive, from Eq. 1, that the energy variation for each star is 
proportional to the angular momentum variation; 

AE = QpAJ. (4) 

Eq. 4 is not new and it was obtained also by Lynden- 
Bell & Kalnajs (1972) and Sellwood & Binney (2002). The 
process of steady angular momentum transfer between the 
spiral density and a star in resonant motion with the per¬ 
turbation is understood in the following way: the loss (gain) 
of angular momentum AJ by a star at the inner Lindblad 
resonance is accompanied by the loss of energy QpAJ; at 
the same time, the change in orbital energy, relative to the 
circular motion, is QAJ, which is less than QpAJ. The con¬ 
nection between the orbital energy E and angular momen¬ 
tum J is shown by Eq. 2, where for a pure circular motion 
it becomes Ec = J‘^f{2r'^) -|- 'I>o(r), thus any variation in 
orbital energy leads to a variation in angular momentum at 
a rate proportional to the stellar angular velocity Q. Thus 
the amount of energy in radial direction AEr acquired by 
the star appears as non-circular motion: 

AEr = {Qp — Q)AJ, (5) 

We can merge Eqs. 4 and 5 into the following: 

AEr +AE+ QAJ = 2AJQp. (6) 

The challenge here is to find AEr, AE, and AJ for 
a real set of stars in the Galaxy. The first consideration 
that we can make is that the axisymmetrical potential and 
the kinetic energy are larger than the perturbation, thus we 
might ignore the term 4>i(r, 6^) from Eq. 2. This is justifiable 
because the residual velocities of the objects in comparison 
with the pure circular ones already carries information about 
the perturbing field, and another important reason is that 
we are not giving any prior information about the spiral 
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arms. Thus we can rewrite Eq. 2 ignoring the perturbative 
term; 

E^l(u^ + ^)+Mr)- (7) 

When we observe a star we just have access to its en¬ 
ergy and angular momentum at one specihc period, in other 
words, we cannot track a star for a few millions of years to 
see how it will change its energy and angular momentum 
over time. In an idealized case, where we know exactly the 
values of AE and AJ (e.g. in a simulation), we could use 
Eq. 4 to recover flp. However, the problem of not knowing 
the exact values of AE and AJ introduces errors in this 
computation. This is why we make use of Eq. 6, which is 
the combination of Eqs. 4 and 5, forcing both equations to 
be satished giving us a better result. We also average the 
values oi y = AEr + AE -\- QAJ and x = 2AJ for all stars 
inside a bin, which improves our results as explained below. 
It is important to emphasize that we bin the mean radius of 
the stars and not the actual radial position. 

To solve the problem of the unknown AE and AJ, we 
assume that the initial energy Eq and angular momentum 
Jo of each star is the energy of a circular orbit at the mean 
radius Vm (Eq. 21 gives the definition of r^): 

= ( 8 ) 

with 


a circular orbit in which the star could have come from (see 
Section 2.1 for a better explanation). Therefore, inside each 
bin we have a number N of stars within r < Xm < r-|-Ar and 
we average AEi, AJi and AEr- for all the N stars within 
the bin. Thus Eq. 6 can be rewritten as: 

AE = 2HpAJ, (13) 

where AE and AJ give us the energies and the angular 
momentum variation in each bin; 

- 1 ^ 

AE= — ^(AEi-b AE,,-bH.AJi), (14) 

i = l 

and 

1 ^ 

i=l 

Thereby, using Eq. 13, we can recover the value of Hp 
without giving any previous input about the perturbation by 
htting a first degree equation y = ax + b. Where, y — AE 
and X = 2AJ, with a and h as free parameters. The slope of 
this £t give us the value of Hp. 

In Section 2.1 we show how we compute the mean ra¬ 
dius and the rotation curve. In Section 2.2 we carry out tests 
based on simulated particles in order to show our method 
works, and to estimate the expected uncertainty on the re¬ 
trieved flp parameter. 


Jo — (^) 

where 14 is the rotation curve supplied by Eq. 18. Now we 
have all the necessary ingredients to compute the variation 
in energy and angular momentum: 

AEi =Ei-Eo. = ^+^-~^^+A<ho(ri-r,„,), (10) 


AJi = Ji - Jo- = ril4 - r„.14(rmj, (11) 

and 


AE,. = 


U?-U, 


2 

max 


( 12 ) 


The sub-index i refers to each star, 14 and Ui are re¬ 
spectively, the observed circular velocity and the radial ve¬ 
locity with respect to the center of the Galaxy, and U-max 
is the maximum radial velocity, that happens when the star 
crosses r-m.- Both velocities 14 and Ui, are not directly ob¬ 
servables, what we observe directly are the proper motion 
and the heliocentric velocity, to make this transformation we 
follow Johnson & Soderblom (1987). To correct the motion 
due to the local standard rest (LSR) we used the values from 
Schonrich et al. (2010); uq = —11.1 km s“^ and vo = 12.24 
km s“^. 

The next step is to bin the mean radius in intervals 
of Ar = 0.2 kpc from 5 up to 14 kpc. The reason for it 
is that stars with similar mean radius are expected to have 
almost the same energy H, which reflect the initial energy of 


2.1 Assumptions of the method and uncertainties 

We critically discuss the main assumptions of our method 
and their impact on the recovered Qp value. The main as¬ 
sumption of our method is that there is a dominant pattern 
speed that lasts, at least, for a few billion years. Notice that 
all existing methods are based on this premise. 

As stated before, in a real sample we don’t have access 
to the true variation in energy and angular momentum of 
each star. To overcome this problem another main assump¬ 
tion made here is that the initial energy Eq and the initial 
angular momentum Jo, come from a circular orbit placed at 
the mean radius (see Fig. 1). In a pure circular orbit all the 
stars at the same radius have the same energy and angular 
momentum as those quantities depend only on the stellar 
position (see Eqs. 8 and 9). We then split our sample in bins 
of mean radius (as stars with similar mean radius will have, 
according to our assumption, similar initial energy Eq and 
angular momentum Jo). The bin width we chose was 0.2 
kpc, to assure that the Eqs. 8 and 9 do not vary too much 
and at the same time we could have enough stars (at least 
more than 5) in each bin. We next discuss the impact radial 
migration would have on this important assumption of our 
method. 

Minchev et al. (2013) showed that in a simulated disk 
similar to the Milky Way migration is a global process, sig- 
nihcantly affecting the entire disk. How would this affect our 
results, in particular the relation Qp = AE/AJI We expect 
that the impact of migration on our method will be small 
for the following reasons. Minchev et al. (2012) showed that, 
due to conservation of vertical and radial actions, migrators 
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Table 1. Parameters of the rotation curve (Eqs. 18 and 19). 


a 

[km s“^] 

/3 

[kpc] 

7 

[kpc] 

S 

[km s“^] 

e 

[kpc] 

V 

[kpc] 

Ajxirc 
[km s“^] 

r^mrc 

[kpc] 

^mrc 

[kpc] 

250 

120 

3.4 

360 

3.1 

0.09 

20 

8.8 

0.8 


arrive at a new radial location with orbital properties very 
similar to the stars which did not migrate. Therefore, the 
AE and A J values at the final migration time are expected 
to be very similar to those of the local non-migrators. An 
exception to this rule would come from stars currently in the 
process of migration. While those could be a large number at 
high redshift due to the strong effect of external perturbers 
(e.g., large infalling satellites), it should be expected that 
at present these ’’migrators in action” do not constitute a 
significant fraction of the stars found in a given radial bin. 

Another possible source of error can result from kine¬ 
matically hot stars, i.e., stars with high eccentricity, for 
which our main assumption that the initial energy Eq and 
the initial angular momentum Jo, come from a circular or¬ 
bit at the mean radius rm, starts being imprecise. This is 
the main reason for using averaged AE and AJ values for 
stars sharing the same bin in mean radius. Because the stars 
with hot kinematics are mostly old, we adopt a) a sample of 
open clusters for which only 8% have ages above 1 Gyr (see 
Fig. 3), b) a stellar sample confined to the Galactic plane, 
mostly dominated by thin disk stars and c) spiral arms have 
the strongest dynamical effect in the disk midplane and thus 
stars with low vertical oscillations are the best tracers for the 
arms. As we discuss in the Results Section, the fact that we 
obtain very similar results from both samples, show the im¬ 
pact of the above mentioned shortcomings (radial migration 
and stars on eccentric orbits) in our method to be minor. 

Finally, the main source of error in our method comes 
from the computation of AE and A J. To illustrate how the 
errors from AE and AJ affect the measurement of flp, we 
propagate the errors using the Eq. 4 and we derived the 
equation bellow: 




1 

AJ 



AE2 „ 
AJ2 


(16) 


which can be seeing in the errors bar for AE in Fig. 2, for 
A J the magnitude of the errors bar are the same. Thus, from 
Eq. 17 with a = /3 = 0.1 and flp = 23 km s“^ kpc“^ we 
have an error of = 3 km s“^ kpc“^. In a more general 
way (JOp = O.ldflj, km s“^ kpc“^ 

In summary, the oldest the population, the less precise 
is our method. On the other hand, a very young popula¬ 
tion (<10 Myrs), which most probably would not have had 
enough time to interact with the spiral arms, because of the 
transformation of dissipational to collisionless dynamics, can 
also lead to uncertain results. An optimal sample would be 
composed of stars with ages between 50 Myrs and a few 
Gyrs. Finally, as discussed before we expect larger errors 
for larger values of AE and AJ. Hence, larger uncertainties 
should be expected for stars migrating from the resonances 
(i.e. inner or outer Lindblad resonance, ILR or OLR), as 
these stars might have increased significantly their energy 
and and angular momentum. The two samples adopted here 
were chosen with the aim to minimize these effects (see Sec¬ 
tion 3). Finally, we notice that a large number of stars per 
mean radius bin improves our determination of flp. 


2.2 Mean radius and the rotation curve 

In order to find the mean radius we adopt a model for the 
axisymmetric galactic potential that reproduces the general 
behavior of the rotation curve of the Galaxy. We use an 
analytical expression to represent the circular velocity as a 
function of Galactic radius, conveniently fitted by exponen¬ 
tial in the form (units are km s“^ and kpc): 


Vc(r) = nexp 




-I- (5 exp 




It simplify the analysis once we assume that the major 
errors come only from AE and AJ. Here ctaj and ctab are 
the errors from these two variables. Now let’s assume that 
the errors follow the relation; (Jab = aAE and ctaj = /3A J. 
This tell us that the errors are proportional to the own vari¬ 
ation of energy and angular momentum, respectively. This 
makes sense because Eq and Jo become less precise when 
the variation of AE and AJ are larger. Thus, for the errors 
following the given definition we can rewrite the equation 
above as: 

A p 

o-fip = +/3^ = flp\/a‘^ + 13^. (17) 

For a single star the parameters, a and /3, can be larger 
than one, which leads to an error greater than 100% in flp. 
However, the errors in a and j3 are much smaller once we 
use a large number of stars averaged on a particular mean 
radius bin. Indeed, for the test sample a and jj are « 0.1, 


with 


J'mrc(^) — Anirc 


1 

2 


r - Rn 


2 


(19) 


where Amrc is the amplitude and Umrc is the half-width of 
the minimum centered at the radius Rmrc. We verified that 
the adopted depth of the minimum in Eq. 19 does not have 
any measurable effect on the value of flp obtained in this 
work. The rotation curve given by the expressions in Eqs. 18 
and 19 is close to that derived by Eich et al. (1989) and is 
also similar to the ones previously used by, e.g., Lepine et 
al. 2008; ALM; Lepine et al. 2011a. The interpretation of a 
similar curve in terms of components of the Galaxy is given 
by Lepine & Leroy (2000). Table 1 gives the values of the 
parameters chosen to reproduce the rotation curve of the 
Milky Way. For the Galactocentric distance of the Sun, we 
adopt Ro = 8.0 kpc. The circular velocity at Ro resultant 
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X (kpc) 


Figure 1. This figure illustrates a typical orbit, in an inertial 
frame of reference, where the star follows elliptical orbits which 
are confined in a circular region limited by the pericenter (rp) and 
the apocenter (ra). The black lines show the stellar orbit, the red 
line is the position of the mean radius Vm- 

from Eq. 18 is Vo = 210 km s“^, for a peculiar velocity of 
the Sun in the direction of Galactic rotation vo = 12.24 km 
s“^ the velocity with respect to the local standard rest is 
Vlsr — 220 km s“^. For more details and also a theoretical 
description about the minimum close to the solar position 
see Barros et al. (2013). 

As we restrict our study to orbits in the galactic plane, 
the axisymmetric potential can be derived directly from the 
rotation curve: 



We integrated the stellar orbits for 2 Gyr under the in¬ 
fluence of the potential 4>o(r), excluding any perturbations. 
This allowed us to find the apocenter (ra) and pericenter 
(rp) radius. The mean radius is then given by; 




2AJ (kpc^Myr ’) 



2AJ (kpc^ Myr’) 


_ra+rp 
~ 2 ■ 


( 21 ) 


Fig. 1 shows an example of an orbit with the apocen¬ 
ter, pericenter and mean radius of a star in the Galactic 
potential. 


Figure 2. This figure shows the linear correlation between 2AJ 
and AS for three different values of flp. The red dots are the 
values computed from Eqs. 14 and 15, as explained in Sec. 2, the 
dashed red lines are the least-squares fit on these points, where 
the slope give us the value of tip (shown on Tab. 2), and the black 
lines are the linear relation for the original value of Op. 


2.3 Application of the method to modeled data 

In order to test our method we integrate the orbits of 500 
test-particles for 2 Gyr under the influence of a perturbing 
potential with a given value of Op. We chose 500 particles 
to match with the number of stars we have available in each 
sample, approximately. With future data, the increase of 
these number and a better distributed in the Galactic plane, 
could improve the results, as we will discuss later in the 
Sec. 3. The potential and parameters that we used for the 
spiral arms are described in Junqueira et al. (2013) (their 
Eq. 6 with the parameter values given on their Tab. 1). 
The axisymmetric potential $0 comes from Eq. 20, with 
the rotation curve 14 given by Eq. 18, where we set 'l?o(100) 


= 0 at r = 100 kpc to find the constant of integration. 
Initially, the test-particles were distributed between 6 and 
12 kpc with random azimuthal positions, and initial circular 
velocities corresponding to the rotation curve 14. The test- 
particles have initial radial velocities U, given by a Gaussian 
shape in each bin of radius and the half-width is the velocity 
dispersion Ur that follows a radial profile given by Eq. 22. 

fTr = Uro exp ^ ^ , (22) 

with To = 8 kpc at the solar radius and a velocity disper¬ 
sion Grg = 6.5 km/s. These values are compatible with the 
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Table 2. The first column shows the input value for flp and the 
second column shows the value of Qp recovered by our method. 
All values in this table have units of km s“^ kpc“^. 


Input Qp 

Recovered Qp 

23 

23.5±0.9 

26 

26.7±1.0 

30 

29.7±0.7 


amplitudes of the perturbation velocities due to the spiral 
waves found in the literature (e.g. Burton 1971, Mishurov 
et al. 1997, Bobylev & Bajkova 2010). They are also similar 
to the velocity dispersion of the youngest Hipparcos stars 
(Aumer & Binney 2009). 

The second step is to recover the value of Qp from this 
synthetic sample that we generated. To do that we selected 
the particles from different snapshots in order to simulate 
different ages, as we have in a real sample. After that, we 
re-integrate the orbits using only the axisymmetric potential 
to find the mean radius of each particle. Then we use the 
method explained before to compute the Eqs. 14 and 15. In 
Fig. 2 we show the results obtained for different values of 
Qp, where the error bars are the RMS (root mean square) of 
each bin. Table 2 summarizes the results; the first column 
is the input values of flp, the second column shows the re¬ 
covered values, that are given by the slope of a least-squares 
fit shown in Fig. 2, with the respective asymptotic standard 
errors. We can see that our method is able to recover the 
value of Qp down to around 4% precision. 

In Sec. 4 we will apply this method to a sample of open 
clusters and stars from Apogee catalog to compute the Qp 
of the MW spiral arms. 


3 THE DATA 

In this Section we describe the two data sets we adopt to 
illustrate our new proposed method for deriving the pat¬ 
tern speed of the spiral structure of the Milky Way. As dis¬ 
cussed in Section 2.3 we are interested in a sample of stars 
and/or clusters dominated by young-intermediate age ob¬ 
jects, mostly close to the Galactic plane and for which we 
have have good distances, radial velocities and proper mo¬ 
tion information. Here we adopt two samples, which have 
complementary advantages and disadvantages. In this way 
we are able to illustrate the robustness of our method. In¬ 
deed, as we will see in the Results section, the two samples 
lead to essentially the same results, despite their different 
azimuthal and age coverage. We now describe each of the 
samples. 

3.1 Open Clusters 

Open clusters play an important role on the study of Galac¬ 
tic dynamics, because they are mostly concentrated in the 
disc plane. Thus, we can use them to find evidences about 
the kinematic and evolution of the MW’s disc. 

The open clusters that we used to measure the pattern 
speed of the spiral arms belong to the “New Gatalog of Op¬ 
tically Visible Open Glusters and Candidates”, published 


by Dias et al. (2002)(DAML02)^, version 3.3. This catalog 
is an update from the previous ones by, Lynga & Palous 
(1987) and Mermilliod et al. (1995), and contains 2140 ob¬ 
jects with measured parameters such as distance (for 74.5% 
of the sample), age (74.5%), proper motion (54.7%) and ra¬ 
dial velocity (24.2%). In this work we used 513 open clusters 
from this catalog, which have distance, proper motion and 
radial velocity available simultaneously. Fig. 3 shows the age 
distribution of our OC sample. It can be seen that most of 
the objects are confined in the 10-1000 Myr age range which 
is an ideal age range for applying our method (see discus¬ 
sion in Section 2.3). Finally, the mean radius distribution of 
our OC sample is shown in Fig. 4 (solid black line). The OC 
sample is concentrated in the 7-10 kpc mean radius range 
where the percentage of young stars coming from the res¬ 
onance regions is expected to be small (see Minchev et al. 
2013, 2014) .For the OCs, Paunzen & Netopil (2006) show 
a limit of 20% in the errors for distances, which are similar 
to APOGEE HQ sample, and the errors in proper motion 
for 90% of our OCs sample are less than 1 mas/yr. The 
uncertainties in radial velocities are less than 5%. 


3.2 APOGEE 

The Apache Point Observatory Galactic Evolution Exper¬ 
iment (APOGEE; Allende Prieto et al. 2008; Majewski & 
The SDSS-IH/APOGEE Collaboration 2014), is one of the 
four Sloan Digital Sky Surveys III (SDSSHI; Eisenstein et 
al. 2011). Recently, Anders et al. (2014) have defined a sub¬ 
sample of the first year of APOGEE data (as part of data 
release 10, DRIO; Ahn et al. 2014). The selection criteria for 
what the latter authors named their APOGEE High Qual¬ 
ity Giant Sample are summarized in Table 1 of Anders et 
al. (2014). 

Starting from a similar sample, we have selected stars 
that stay on the galactic plane (i.e. with their maximum 
vertical orbital amplitude, z^ax, below 0.2 kpc. Moreover, 
we required a combined proper-motion error below 4 mas/yr 
< 4 mas/yr and < 4 mas/yr, mean proper mo¬ 

tion error in right ascension and declination, respectively). 
The final sample resulted in 559 stars from DRIO which is a 
subsample from what Anders et al. (2014) named their gold 
sample. The distances were computed using the distance 
code of Santiago (priv. communication - see also Santiago 
et al. 2015). The mean uncertainties of distances and proper 
motions for the APOGEE DRIO sample are shown by An¬ 
ders et al. (2014), where their gold sample have a threshold 
in uncertainties of 20% in distances and 4 mas/yr in proper 
motion. Also for the APOGEE sample, the uncertainties in 
radial velocities are less than 5%. Their mean radius dis¬ 
tribution are shown in Fig. 4 (dashed red histogram), and 
turned out to be similar to that of our sample of OCs. The 
main difference is that the APOGEE giants span, most prob¬ 
ably a larger age range. Indeed, we expect stars with mean 
radius in the 7-10 kpc range to be predominantly of ages 
between 1 and 6 Gyrs (see Minchev et al. 2014). For consis¬ 
tency, for the final sample of APOGEE stars we recompute 
the mean radius r„j using the rotation curve given by Eq. 18 


^ http://www.astro.iag.usp.br/ocdb/ 
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Figure 3. Logarithmic distribution of the ages for the OC’s sam¬ 
ple. 




(tpc) 


Figure 4. This figure shows the mean radius distribution of the 
OC sample (solid black curve) and red giant branch stars from 
APOGEE (dashed red curve). Both samples have similar distri¬ 
bution, as well as almost the same number of objects per bin. 


instead of the one from Anders et al. (2014). However, in 
both calculations the mean radii are very similar. 

Finally, in Fig. 5 we show the spatial distribution of 
both samples on the X-Y Galactic plane. The first thing 
that can be noticed is that the OCs are more homogeneously 
distributed in the azimuthal direction than the APOGEE 
sample adopted here. The second difference is that we have 
more OGs in the inner part of the disk, while the APOGEE 
sample is more concentrated in the outer parts. This will 
certainly be improved once APOGEE-2 data will be avail¬ 
able. As we will see in the next Section, despite these main 
differences, the pattern speed computed with both samples 
turned out to be very similar. 


4 RESULTS AND DISCUSSION 

In Sec. 2 we described and tested a new method, without 
invoking any prior information about the spiral arms, which 
proved to be useful to constrain the value of the pattern 
speed within an error of ~ 1 km s“^ kpc“^. Here we show 
the results that we obtain by applying our new method to a 
sample of disk stars and open clusters, described in Sec. 3. 


Figure 5. This figure shows the distribution in the X-Y galac¬ 
tic plane, with the Sun (yellow dot) locate in (0,0), for our 
adopted OC sample (black crosses) and red giant branch stars 
from APOGEE (red open circles). The Galactic center (GC) is 
pointed toward negative values in the Y-direction indicated by a 
black dot. We can see that the OC sample has a more wide az¬ 
imuthal distribution than the giant branch stars from APOGEE 
(see text for more details). 


4.1 The value of flp and the corotation radius 

Our main results are summarized in Fig. 6, where we plot 
2AJ versus AE. The slope of the observed linear correla¬ 
tion gives the value of Op (shown in Tab. 4), as described 
by Eq. 13. The upper, middle and bottom panels show the 
results we obtain using the OC, red giants and the combi¬ 
nation of the two, respectively. The color gradient indicates 
the center of the bin in radius in which each point. The value 
that we find by combining the two samples is Op = 23.0±0.5 
km s“^ kpc“^. For Op = 23 km s“^ kpc“^ the corotation 
radius in our model is situated at Ter ~ (1.09 ± 0.04)Ro 
which is in agreement with results found by Dias & Lepine 
(2005) with Ter = (1.06±0.08)Ro. The fact that the Galaxy 
presents a well defined corotation radius supports the idea of 
a dominant pattern speed, at least, for an old stellar popu¬ 
lation (see further discussion in Section 4.3). It’s interesting 
to notice that the values obtained for Dp from the fit, in 
Fig. 2 or Fig. 6, have unit of Myrs“^. However the usual 
unit given for the pattern speed is in km s“^ kpc“^, thus to 
transform from (Myrs”'^) to (km s“^ kpc“^) we must divide 
by a factor of 0.00102. 

We should keep in mind that our quoted uncertainties 
take into account only the fitting procedure. Systematic and 
intrinsic errors (e.g. assumptions made on the method and 
variation of the rotation curve because of uncertainties on 
the solar parameters Ro and Vo) were not taken into ac¬ 
count. However, we can estimate that by varying the values 
of Ro and Vo between 7 — 8.5 kpc and 199 — 245 km s“^ re¬ 
spectively. These variations would propagate into a ~2 km 
s“^ kpc“^ change on the final value of Dp. Piffl et al. (2014) 
recently found Vo = 240 km s“^ for Ro = 8.3 kpc, which is 
compatible with our adopted values and hence will change 
our estimate of Dp within the intrinsic errors. In Sec. 2.1 
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Table 3. Compilation of some estimations of Qp found in the literature. 


flp (km s ^ kpc 

Method“" 

Objects 

Reference 

19.1±3.6 

1 

Cepheids 

Mishurov et al. (1979) 

30±0.7 

1 

O and B-type stars and Cepheid 

Fernandez et al. (2001) 

20±2 

2 

Open clusters 

Amaral & Lepine (1997) 

24±1 

2 

Open clusters 

Dias & Lepine (2005) 

20.3±0.5 

2 

Sample of runaways and early-type stars (from Hipparcos) 

Silva &; Napiwotzki (2013) 

21.2±1.1 

2 

Only early-type stars 

Silva & Napiwotzki (2013) 

18.1±0.8 

3 

Hipparcos subsample 

Quillen & Minchev (2005) 

18.6±0.3 

4 

RAVE survey 

Siebert et al. (2012) 


“"method 1 = kinematic model, method 2 = Birthplace technique, method 3 = Orbital analysis of moving group in the (u-v) plane, 
“"method 4 = Spiral perturbation to reproduce the gradient in the mean galactocentric radial velocity. 


we discussed more about the intrinsic uncertainties of our 
method. 

As we can see in the last panel in Fig. 6 the combination 
of both data, APOGEE + OCs, has a better distribution 
in 2AJ, this is due to a better spatial distribution when 
we combine the data. It happens because the stars have 
different values of radial velocities (and they do not belong to 
a group in a U-V velocity diagram) given us a better average 
of the angular momentum variation in an annular region of 
the disk. Thus, we expected that with the new coming data, 
from APOGEE-2 (now starting with SDSSIV), our results 
will improve, since it will increase the number of stars and 
in a way that the distribution in azimuth will become more 
homogenous. Another thing we can notice in this figure is 
that the points with higher |2AJ| and \AE\ come from mean 
radius around 6 or 12 kpc, that are very close to the IRL 
and OLR resonances and as we discussed before this regions 
warm up the disk. 


4.2 Comparison with previous reported results in 
the literature 

One of the most important parameters in studying the spiral 
structure is its pattern speed Qp. Although the fundamental 
nature of the spiral arms is not fully understood it plays an 
important role in galactic dynamics (e.g. Antoja et al. 2009; 
Lepine et al. 2011b; Quillen et al. 2011; Minchev et al. 2012; 
Sellwood 2014) and its pattern speed is a fundamental pa¬ 
rameter that drives all the resonances in the disc. Table 3 
summarizes the several previous attempts made in the lit¬ 
erature to estimate the value of flp. Gerhard (2011) made a 
review from the values found for the pattern speeds in the 
Milky Way and he end up with a range between flp ~ 17—30 
km s“^ kpc“^. 

As it is clearly seen from Table 3 and 4, our results 
are in agreement with studies that suggest a pattern speed 
between 20-25 km s“^ kpc“^. However, Quillen & Minchev 
(2005) and Siebert et al. (2012) found values bellow 20 km 
s“^ kpc“^, which are not in agreement with our results even 
with error bars around ~ 3 km s“^ kpc“^. In the literature 
we see that higher values for Hp are preferred by open clus¬ 
ter birthplaces while hydrodynamical simulations and phase 
space substructures favor slower pattern speeds. Thus, since 
Quillen & Minchev (2005) and Siebert et al. (2012) anal¬ 
ysis are based on moving groups and the velocity gradient 
both close to the solar neighborhood, which are substruc- 


Table 4. This table contains our results for the pattern speed Qp 
obtained from two different samples and the combination of all 
them. 


$lp(km s ^ kpc ^) 

Sample 

24.0±1.0 

Open clusters 

23.3±0.6 

Apogee 

23.0±0.5 

Apogee -h Open clusters 


tures and can be associated only with one spiral structure 
(as e.g., Perseus arm or even a local spur arm), it could ex¬ 
plain why they found lower values for Q,p. The discrepancies 
of values found in the literature are discussed in Sec. 4.3, as 
a possible contamination by multiple spiral patterns, which 
are difficult to be taken into account and can lead to a sys¬ 
tematic errors that explain the wide range of values found 
for Qp, depending on the tracers and the methods that were 
used to estimate it. 

The value that we found for the dominant MW spiral 
pattern speed is also in agreement with the values found 
in many external galaxies. For example, Scarano & Lepine 
(2013) found flp to have a distribution concentrated around 
24 km s“^ kpc“^. However, the fundamental nature of the 
pattern speed is still not clear, which requires more theoret¬ 
ical work to be fully understood. 

4.3 Possible contamination by multiple spiral 
patterns 

Vallee (2014) uses many tracers to probe the spiral struc¬ 
ture of our Galaxy and concludes that it has a four-armed 
spiral pattern. However, only two of these may be present in 
the density distribution of old stars (see, e.g. Drimmel 2000; 
Martos et al. 2004). Multiple spiral patterns could possibly 
be a source of errors when we try to determine only one 
pattern speed. Naoz & Shaviv (2007) measured the pattern 
speed for: Sagittarius-Carina and they found a superposi¬ 
tion of two pattern speeds with Q,p — 16.5 ± 2.0 km s“^ 
kpc“^ and Q,p = 29.8 ± 2.0 km s“^ kpc“^, Perseus arm 
Hp = 20.0 ± 2.0 km s“^ kpc“^ and Orion Hp = 28.9 ± 2.0 
km s“^ kpc“^. Some models also support multiple pattern 
speeds in order to explain radial migration in discs of galax¬ 
ies (Minchev & Quillen 2006; Minchev et al. 2012; Grand 
et al. 2014). However others studies suggest that the MW 
has a corotation radius well established, situated close to 
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3.007 

-0.18 -0.14 


-0.1 -0.06 -0.02 
2Aj (kpc^Myf') 


0.02 0.06 


The fact that the Galaxy presents a well defined corotation 
radius supports the idea of a dominant pattern speed, at 
least, for an old stellar population. Therefore, one way to 
avoid a possible contamination by different pattern speeds 
is to split the sample into old and young stars. 

In the future we will apply our method to a N-body 
simulation sample with multiple pattern speeds to check if 
we are able to distinguish multiple arms and/or analyze the 
influence of small arms on the dominant ones. 


5 CONCLUSIONS 
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Figure 6. This figure shows the linear correlation between 2AJ 
and AE for a real sample. The dots are the values computed 
from Eqs. 14 and 15, as explained in Sec. 2, and the color code 
represent the bin radius in which each dot was computed. The 
black lines are the least-squares fit on these points, where the 
slope gives us the value of tip. The upper panel shows the result 
for the OC sample; the middle shows the results for the APOGEE 
sample and the last one is the combination of both, APOGEE -|- 
OCs samples. The values of Qp found for each sample are shown 
on Tab. 4. 

the solar radius (Marochnik et al. 1972; Creze & Mennessier 
1973; Mishurov & Zenina 1999; Dias & Lepine 2005; Amores 
et al. 2009, among others). It would be possible only if we 
have a dominant pattern speed rotating rigidly, neverthe¬ 
less this do not exclude small structures to exist, that may 
rotate with different angular velocities. For Qp = 23 km 
s“^ kpc“^ the corotation radius in our model is situated at 
rcr = (1.09 ± 0.04)Ro which is in agreement with results 
found by Dias & Lepine (2005) with r^r = (1.06 ± 0.08)Ro. 


In this work, we proposed a new method to derive the spiral 
pattern speed of the MW based on the interaction between 
the spiral arms and the stellar objects. In this method we 
do not need any prior information about the spiral arms, as 
for example its shape and location. In addition, we do not 
need to be concerned with the stellar ages, which allows the 
use of data from large Galactic surveys. 

The assumption we make in our method is that the 
initial energy and angular momentum of the objects can be 
approximated as the circular orbit at the mean radius. This 
approach introduces a natural error on the order of 1~2 km 
s“^ kpc“^ which is equivalent, or even smaller, than the 
available methods to constrain the value of Up. 

Using a sample of open clusters and red giant stars 
from the APOGEE DRIO (Ahn et al. 2014) we have found 
Qp = 23.0 ± 0.5 km s“^ kpc“^, which is compatible with 
other values in literature and placed the corotation radius 
at Ter = 8.74 kpc for a solar position Rq = 8 kpc. We have 
to stress again that the given error for Dp here is just due to 
the RMS from the fitting procedure. A more realistic error 
estimate which also takes into account the errors intrinsic 
to our method should be around 2 km s“^ kpc“^ (intrin¬ 
sic method error -|- RMS). Systematic errors and errors due 
to other sources of perturbation (as discussed in Sec. 4.3) 
are even more difficult to estimate, and could most proba¬ 
bly explain the range of values for Dp, between 17 — 30 km 
s“^ kpc“^, found in the literature. Further studies, using 
N-body simulations data, are needed to check the effective 
influence of multiple spiral arms on the determination of Dp, 
assuming a constant pattern speed. 

The new method for estimating the spiral pattern speed 
presented here can be tested with the large amounts of cur¬ 
rently available data of ever increasing quality from large 
Galactic surveys, such as RAVE (Steinmetz et al. 2006), 
SEGUE (Yanny et al. 2009), APOGEE (Allende Prieto et 
al. 2008; Majewski & The SDSS-III/APOGEE Gollabora- 
tion 2014), GES (Gilmore et al. 2012), and in the near fu¬ 
ture - Gaia (de Bruijne 2012), 4MOST (de Jong et al. 2012) 
and WEAVE(Dalton et al. 2012). 
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